APPENDIX 



.MODULE/ABS=0LX1; 
.pagelength 44; 

r This software is written in the native assembly language 
of the Analog Devices 2181 Digital Signal Processor. 
The purpose of this code is to oversample bead flu rescence 
in each of 5 consecutive channels t isolate single bead events 
and report the fluor area in each channel. The program 
double buffers the data for the 167 micro, and uses the 
following interrupts for interprocessor signalling 
From 167 to IRQ this program: 
p2.0_IRQ0 go to idle state 

P2.1 IRQ1 new data in control area 

P 2 2 — Fl begin capture (this is a state, not IRQ) 

p2.3_PWRDN nmi begin execution at PM 0x2c 

From this program to IRQ 167: 
FL0_EX0IN buffer 0 is ready 

FL1_EX1!N buffer 1 is ready 

FL2 EX2IN error occurred during capture state 

All PM and DM are accessible from the 167 thru IDMA so 

mem areas are asigned and restricted. 

PM is loaded from EPROM but may be rewritten bv 

1) IRQ0(p2.0) 

2) rewrite with vector to first instruction at 0x2c 

3) PWRDN nmi (p2.3) 

DM allocation (in hex) is: 

0-ff reserved for 2 1 8 1 

100- Iff DSP status area, read only to 167 

200-2ff DSP control area, read all/ write restricted to 167 

300-3ff DSP command area, read/write 167, read / restricted write to DSP 

400-9ff buffer 0, 1536 words, 384 max block size 

aOO-fff buffer 1 

1 000-3dff reserved for DSP raw A/D data 



7 

.CONST PFLAGS=0X3FE5; 
r DSPs variable area V 

.VAR/dm/ABS = 0x90 inbkgnd[10]; 
.VAR/dm/ABS = OxaO x8bkgnd[10]; 
.VAR/dm/ABS = OxbO histbk[10]; 
r DSP status */ 

52 7, 0X100 miSSCnt ' M ' buf0cnt ' »2. buflcnt otrcnt,xx4 I afreg,xx5,nobgsub; 

fmrn^SS!^ 0X200 ^ X31 ' bUf1SiZ ' X32 ' thr6Sh ' X33 ' mCre ^ ™ 4 ' dock ' 
r DSP command V 

.VAR/dm/ABS = 0X300 address,xxb1 t hiword,xxb2,lobyte f xxb3,msq- 
.VAR/dm/ABS = 0X400 outbufO; 9f 

.VAR/dm/ABS = OXaOO outbufl ; 
r raw data area V 

.VAR/dm/abs = 0x1000 hist[0x1000] ; 
.VAR/dm/circ/abs = 0x2000 workbuf[0x1e00] ; 



r interrupt vectors 7 



dis ints; 
jump begin; 
nop; 
nop; 

jump irqf); 
nop; 
nop; 
nop; 

rti; 
nop; 
nop; 
nop; 

rti; 
nop; 
nop; 
nop; 

rti; 
nop; 
nop; 
nop; 

rti; 
nop; 
nop; 
nop; 

rti; 
nop; 
nop; 
nop; 

jump bdmadone; 

nop; 

nop; 

nop; 

jump procnewctl; 

nop; 

nop; 

nop; 

jump doidle; 
nop; 
nop; 
nop; 

jump tim_exp; 
nop; 
nop; 
nop; 



Treset irqV 

Tread 4 fluorescence measurements*/ 
/*irqL1 7 



/*irqLO 7 



/*sp0tx7 



TspO rxV 



Hrq E 7 



/*irq BdmA7 



/*irq1-new data in control area7 



rgo to idle mode, wait for nmi7 



/*timer7 



softload: 



m0=0; 
10=0; 

iO= A address; 



ax0=-1; 

dm(iO,mO)=axO; 

11=0; 

12=0; 

15=0; 

ax0=0xd000; /* set up restart params 7 
dm(0x3fef)=ax0; 
i1= A hiword; 
i2= A lobyte; 
Idpgm: ay0=dm(i0,m0); 
ar=pass ayO; 

if It jump Idpgm; /* negative number means no new program yet 7 
i5=ay0; r point to addr of new instruction 7 

ax0=dm(i2,m0); r get low byte 7 
px=ax0; 

ax0=dm(i1 ,m0); r get hi word 7 
pm(i5,m5)=ax0; Tm reg doesnt matter 7 
ax0=-1; 

dm(iO,mO)=axO; 

ar=pass ayO; A* if 0 then its the last instruction 7 

ifeqjumpldend; 

set M2; 

nop; 

reset M2; 

jump Idpgm; 
Idend: idle ; 

nop; 

nop; 

nop; 

nop; 

nop; 

nop; 

nop; 

nop; 

nop; 

nop; 

nop; 
nop; 
nop; 
nop; 
nop; 
nop; 
nop; 
nop; 
nop; 
nop; 
nop; 
nop; 
nop; 

/* initalization sequence 7 

begin: imask=0; 

icntl=0x7; /*no nesting, irq2, 1 ,0=edge7 

dis timer; 

10=0; 

11=0; 

12=0; 

13=0; 

17=0; 

16=0x1 eOO; 

mstat=0x1 0; ^integer placement7 

ax0=0; 



» 



dm(lena)=axO; 

m0=2; 

m1=-1; 

m2=1; 

m3=0; 

m4=-20; 

m5=1; 

m6=5; 

m7=0; 

M=0x3fff; 

dm(i1,m1)=0; 

dm(i1 ,m1)=0x3ff8; r atod is ZERO wait state 7 

dm(i1,m1)=0; HperiodV 

dm(i1 t m1)=0; /*tcount7 

dm(i1,m1)=15; rtscale to 1 us per cntV 

dm(i1,m1)=0; /*sp07 

dm(i1,m1)=0; 

dm(i1,m1)=0; 

dm(M.m1)=0; 

dm(i1,m1)=0; /*3ff6*/ 
dm(i1,m1)=0; 
dm(i1,m1)=0; 
dm(i1,m1)=0; 

dm(i1 ,m1)=0x4000; rena internal clock genV 

dm(i1,m1)=15; /*ser clock out to 1.024 MHzV 

dm(i1,m1)=0; /*addr3ff07 

dm(i1 t m1)=0xd000; r set up power down options 7 

M=0x3fe5; 

dm(i1 ,m2)=0x0; Tset pfO-7 to 0 7 
dm(i1 f m1)=0x78ff; TpfO-7 are outputs*/ 
reset AO; 
reset 111; 
reset f!2; 

reset flag_out; /* set flags for reg control 7 

axu=0; 

dm(msg)=ax0; 

dm(mcready)=axO; 

dm(bcntl)=axO; 

dm(bcnth)=axO; 

dm(nobgsub)=axO; 



cntr=0xc00; 
iO= A outbufD; 
do clrmem until ce; 
clrmem: dm(iO,m2)=axO; 
14=0x1 eOO; 



I* let system stabilize 7 

cntr=1000; 

do bpirq until ce; 

cntr= 16000; 

do bpirq 1 until ce; 
bpirq 1: nop; 
bpirq: nop; 



irqfl: 



my1=i (0); /* init for first readapds 7 
i1 = 0x3f 5; 

dm(i1 ,m2)=0x0; Tset pfO-7 to 0 */ 
call lightsoff; 

ifc=0xff; Tclear all pending irqsV 

jump mainbody; r for testing only V 

jump waitmc; /* wait for micro to init data areas 7 

r INTERRUPT SERVICE ROUTINES- 7 



my1=io(1); r read 0 setup 1 7 
dm(i4,m5)=my1; 

my1=io(2); r read 1 setup 2 7 
dm(i4,m5)=my1; 

my 1 =io(3); r read 2 setup 3 7 

r f()r new boafd change tQ iQ(4) be|QW v 

r for Q|d bQard change tQ jQ(3) be|Qw #/ 

dm(i4,m5)=my1; 

my1 =io(4); /* read 3 setup 4 V 
dm(i4,m5)=my1 ; 

my 1 =io(0); r read 4 setup 0 V 
dm(i4,m5)=my1 f af=af+1; 



rti; 



procnewctl: 
rti; 

doidle: 

rti; 

bdmadone: 
rti; 



tim_exp: 

rti; 



This loop is executed while NOT capturing. It only responds to 2 interrupts 

IRQ 1 means that there is new data in the control area so reinit all variables 

IRQ 0 means that there is new code in dm that must be moved to pm and executed 



V 



Mainbody: 

if flag jn call capture; 

axO=dm(iena); /*tum on lasers V 

ar=setbit 7 of axO; Tenable overflow bit of a/d converter V 

dm(pflags)=ar; 



imask = 0x206; fallow idle or n wctlmsgV 

enaints; 

nop; 

nop; 

nop; 

nop; 

iO=*workbuf; 

ax0=dm(i0,m2); r write background in same order as recbuf V 
dm(bkgnd)=axO; 

axO=dm(iO,m2); 
dm(bkgnd+2)=ax0; 

axO=dm(iO,m2); 
dm(bkgnd+4)=ax0; 

axO=dm(iO,m2); 
dm(bkgnd+6)=ax0; 

axO=dm(iO,m2); 
dm(bkgnd+8)=ax0; 



jump mainbody; 



I* All subroutines here at the end 7 



capture: 

14=0x1 eOO; 
15=0x1 eOO; 

axO=dm(lena); /*tum on lasers */ 

ar=clrbit 7 of axO; r disable overflow bit of a/d convt 7 

dm(pflags)=ar; 

dis ints; 

call clitiist; 

af=pass 0; 

»4= A workbuf; 

i3= A hist; 

i5= A workbuf; 

i7= A x8bkgnd; 

imask=0x200; 

ena ints; 

cntr=4; r do four background histograms */ 
do bigloop until ce; 

cntr=4000; r take 4000 samples for histogram V 
do bldhist until ce; 
pace: ar=pass af; /* any samples to proc V 
if le jump pace; 

ax0=dm(i5,m6); r get samp V 
ar=ax0-0x1 000; r range check V 
if ge jump bldhist; 
ayO= A hist; r base addr */ 



bldhist: 



ar=axO+ayO; 
i3=ar; 

ay1=dm(i3,m3); 
ar=ay1+1; 
dm(i3,m3)=ar; 
nop; 



r base + offset 7 
/•point into hist 7 
/* get current val 7 



r now find most frequent occurrance in hist 7 



i3= A hist; 

cntr=Oxfff; 

ax0=dm(i3,m2); 

dm(svi3)=i3; 

ayO=axO; 



/* get first entry */ 

r ayO holds highest val so far 7 



do findmax until ce; 
ax0=dm(i3,m2); 
ar=axO-ayO; 
if le jump findmax; 
ayO=axO; 
dm(svi3)=i3; 
findmax: nop; 
ax0=dm(svi3); 
ayO= A hist + 1; 
ar=axO-ayO; 
si=ar; 
ayO= A hist; 
ar=ar+ayO; 



/•new high val found 7 
r save offset +1 (post inc) 7 

/*getptrtobest+1 7 
r since i3 is plus 1 */ 
Tcompensate for base V 
P save bkgnd for shift later */ 



i3=ar; 

ay0=dm(i3 t m1); 

ax0=dm(i3,m0); 

ay1=dm(i3,m1); 

ax1=1; 

ar=ax0-ay1; 

if It jump poside; 

ax1=-1; 

ay1=ax0; 



Tcalc index to peak */ 

Tpoint to peak in hist V 
/*get pk, pnt to previous •/ 
/*get prev, pnt to other side of pk */ 
r pnt back to peak */ 

r which pk is second highest V 

r int pk will be less V 



9os\6e: r now ay1 has val of 2nd peak V 

sr=lshift si by 3 (lo); r mult by 8 */ 
si=ayO; 

ayO=srO; TayO is 8x bkgnd and is ready for adj */ 

sr=lshift si by -2(lo); r get 1/4 of peak height */ 
mtlp: ar=ay1-srt); r reduce 2nd pk by 1/4 of main 7 
if It jump intend; 
ay1=ar; 

ar=ax1 +ayO; Tadjust background + or - 7 

ayO=ar; 

jump intlp; 

intend: dm(i7,m5)=ay0; A* save final bkgnd 7 
call clrhist; 

bigloop: modify(i5,m5); Tpoint to next channel 7 

i3= A x8bkgnd; 

i1= A bkgnd; 

cntr=4; 

do div8 until ce; 

si=dm(i3,m2); /*get avg background , inc +27 
sr=lshift si by -3 (lo); 
div8: dm(i1,m0)=srf); 



af=pass 0; 
axO=0; 

dm(misscnt)=axO; 

dm(buf0cnt)=ax0; 

dm(buf1cnt)=ax0; 

dm(otrcnt)=axO; 

dm(actbuf)=axO; 

axO=dm(clock); 

dm(0x3ff1)=ax0; r set sample elk freq 7 

i4= A workbuf; 

i5 =A workbuf; 

ifc=0xff; 

imask=0x200; I* ena fl ints 7 
i0= A outbuf0; 
resetc: If not flag Jn jump capexit; 

ax0=0; 

ar=pass af; 
dm(afreg)=ar; 
ena ints; 

newsamp: if not flag_in jump capexit; 

ar=pass af; r load axO with # of new samps */ 

ar=ar-64; r are we 64 records ahead*/ 

if le jump newsamp; 

af=af-1 ,ay0=dm(i5,m6); /*get ssc, inc to next record 7 
axO=dm(thresh); 

ar=ax0-ay0; V above thresh then its a bead */ 
if ge jump newsamp; 

/* got a bead V 

dm(pki5)=i5; 
dm(fstart)=i5; 

ax1=ay0; r init peak ssc V 

ay 1=0; r init width count */ 



moresamp: 

af=af-1 i ay0=dm(i5 i m6); /* next next ssc, dec rec count */ 

ar=ay 1 + 1 ; r inc width count 7 

ay 1 =ar,ar=ax1 -ayO ; r is new ssc higher */ 

if ge jump ckdone; 

dm(pki5)=i5; 

ax1=ay0; r save peak ssc V 

ckdone: ar=pass af; r out of records V 

if eq jump bubble; r if 0 then bubble 7 

ar=axO-ayO; r back below thresh yet? 7 
if le jump moresamp; 
i6=dm(pki5); 

modify(i6,m4); r since post inc, sub 20 7 
dm(wid)=ay1 ; 

r build a record with peaks or estimates of peaks 7 

I2=0x1e00; 

i2=i6; 

call i2sum; 

dm(recbuf)=ay1 ; /*store ssc in recbuf*/ 

modify(i6,m5); r add 1 to point to fl2 pk -1 5 7 

i2=i6; 

call i2sum; 

dm(recbuf+1 )=ay1 ; /* store c!1 sum 7 



modify(i6,m5); r add 1 to point t A3 pk -1 5 7 

i2=i6; 

call i2sum; 

dm(recbuf+2)=ay1 ; r store cl2 sum 7 

modify(i6,m5); r add 1 to point to U pk -1 5 7 

i2=i6; 

call i2sum; 

dm(recbuf+3)=ay1 ; r store cl3 sum 7 

modify(i6 f m5); r add 1 to point to f!1 pk -1 5 7 

dm(svi6)=i6; Malign point with red laser 7 

modify(i6 ( m4); r back up 20 more to fl1 peak -35 7 

modify(i6,m4); /* back up 20 more to fI1 peak -55 7 

/*modify(i6,m4); back up 20 more to f!1 peak -75 7 

r since fi1 is a different laser, refind peak to allow for alignment 
differences, search from -15 to +16 7 

ax0=dm(i6,m6); r get 1st fl1 sample 7 
m7=160; 

m7=350; rfor new filter at 50 khz 7 

modify(i6,m7); 

m7=0; 

ay0=dm(i6,m6); 

an=ax0+ay0; 

sr=lshift ar by -2 (LO); 

ay0=sr0; 

ar=ar+ay0; 

dm(x8bkgnd+4)=ar; 



i6=dm(svi6); r test using no peak find 7 
m7=40; 

modify(i6,m7); Tmove right 8 samples due to skew of 50k filter 7 
m7=0; 

cntr=7; 

modify(i6,m6); 
dm(pki5)=i6; 

ay0=dm(i6,m6); r get 1st f!1 sample 7 

a*1 =ay0; r move to old fl1 val 7 

do rplpk until ce; 

ay0=dm{i6,m6); 

ar=ax1 -ayO; r is it a new peak?7 
if gt jump rplpk; 

dm(pki5)=i6; I* store addr +5 of new peak (post inc) 7 
ax1=ay0; 
rplpk: nop; 

i6=dm(pki5); 

modify(i6,m4); r sub 20 to back up 3 samples from pk (was +5)7 
modify(i6,m4); /* sub 20 to back up 4 samples 50 k version7 
i2=i6; 

call i2sum14; 

sr=lshift si by -3 (HI); r shift carry 7 

sr=sr or Ishift ar by -3 (LO); 

dm(recbuf+4)=sr0; Tstore rp1 7 

'2=0; r restore for non circ buffr use 7 

r ok, we have singlet and recbuf is built 7 

i7= A recbuf; 

i3= A x8bkgnd; 

cntr=4; 



do subback until ce; 
ax0=dm(i7,m7); r get val inc 0 7 
ay0=dm(i3,m2); 
an=axO-ayO; 
if It ar=pass 0; 
subback: dm(i7,m5)=ar; 
axO=dm(nobgsub); 
ar=pass axO; 
if ne jump recbdone; 
axO=dm(i7,m7); r get val inc 0 7 
ay0=dm(i3,m2); 
ar=axO-ayO; 
if It ar=pass 0; 
dm(i7,m5)=ar; 

/**""~Sub rfl bg from rp1 ****—-/ 



recbdone: 

r record bead count 7 
ay1=dm(bcntl); 
ar=ay1+1; 
dm(bcntl)=ar; 
ar=dm(bcnth); 
ar=ar+C; 
dm(bcnth)=ar; 



i2= A recbuf; 
ay1=dm(actbuf); 
none=pass ay1; 
if ne jump usebufl; 



ay1=dm(buf0cnt); 
none=pass ay1; 
if gt jump noiOinit; 
iO= A outbufD; 

noiOinit: 

ax1=dm(buf0siz); 

none = ax1 -ay1; 

if ne jump noover; 

ay1=dm(misscnt): 

ar=ay1+1; 

dm(misscnt)=an 

jumpresetc; 
noover: cntr=5; 

do movoutO until ce; 

ar=dm(i2,m2); 
movoutO: dm(iO,mO)=ar; 



r is this first of block 7 
r check for overflow 7 

r store ssc inc by 2 for even addr7 



ayO=dm(buf0cnt); 

ar=ay0+1; 

dm(bufOcnt)=ar; 

ay0=dm(buf0siz); r enough records to irq 167? 7 

ar=ar-ayO; 

if ne jump resetc; 

set AO; 

n °P; /* c167 needs 50 ns min to get edge 7 

nop; 

reset AO; 

ay1=1; 

dm(actbuf)=ay1; r change to buffer 1 7 



jump resetc; 



usebuff: 

ay1=dm(buf1cnt); 
none=pass ay1; 

if gt jump noil init; P is this first of block 7 

iO= A outbuf1; 

noil init: 

ax1 =dm(buf 1 siz); f check for overflow */ 

none = ax1 - ay1; 

if nejump nooverl; 

ay1=dm(misscnt); 

ar=ay1+1; 

dm(misscnt)=ar; 

jump capture; I* major overflow reset and restart capture 7 
nooverl : cntr=5; 
do movoutl until ce; 

ar=dm(i2,m2); 

movoutl : dm(iO,mO)=ar; /* store ssc inc by 2 for even addr7 

ay0=dm(buf1cnt); 

ar=ay0+1; 

dm(buf1cnt)=ar; 

ay0=dm(buf 1 siz); /* enough records to irq 167? 7 

ar=ar-ayO; 

if ne jump resetc; 

set H1; 

nop; r c167 needs 50 ns min to get edge 7 

nop; 

reset fI1; 
ay 1=0; 

dm(actbuf)=ay1 ; I* change back to buffer 0 7 
jump resetc; 

capexit: 

i1 = 0x3fe5; 

dm(i1 ,m2)=0x0; Tset pfO-7 to 0 7 

ifc=0xff; 
rts; 



i2sum14: 

m0=5; 

ar=pass O,ay1=dm02,mO); 
si=0; T use for carry 7 
dm(svcry)=si; 
cntr=20; 

do grate4 until ce; 
ar=ar+ay1 ,ay1=dm(i2,m0); 
if not acjump grate4; 
dis ints; 
ena sec_reg; 
ay0=dm(svcry); 
ar=ay0+1; 
dm(svcry)=ar; 
dis sec_reg; 
ena ints; 
grate4: nop; 

ay1=ar; 
si=dm(svcry); 
m0=2; 
rts; 



i2sum: mO=5; 

ar=pass 0,ay1=dm(i2,m0); 
cntr=8; 

do grate until ce; 
grate: ar=ar+ay1 t ay1=dm(i2,m0); 
ay1=ar; 
m0=2; 
rts; 

bubble: 

if not flag jn jump capexrt; 

ar=pass af; r load axO with # of new samps 7 

ar=ar-256; r are we 64 records *4 ahead 7 
if lejump bubble; 

af=af-1,ay0=dm(i5 f m6); A*get ssc, inc to next record 7 
axO=dm(thresh); 

ar=axO-ayO; /* are we below thresh yet 7 
if le jump bubble; r no, keep waiting 7 

jump newsamp; 



lightsoff: 



clrhist: 



zit: 



ax0=0x24d; 
ax0=0x63; /** 
dm(version)=axO; 
rts; 



i3= A hist; 

axO=0; 

cntr=0x1000; 

do zit until ce; 

dm(i3,m2)=ax0; 

rts; 



*test 



********** # 



waitmc: 



toggle fo; 

axO=dm(mcready); 
a repass axO; 
if eq jump waitmc; 
jump mainbody; 



r wait until mc writes non zero to mcready 7 



-ENDMOD; 



